library(Hmisc)
AthleticsData <- read.csv("http://www.statpower.net/R2101/AthleticsData.csv")
attach(AthleticsData)
AD.R <- cor(AthleticsData)
source("http://www.statpower.net/Content/312/R Stuff/AdvancedFactorFunctions.R")
pureCFA.fit <- QuickCFA(AD.R,3,1000,"pureCFA",c("Hand-Eye","Endurance","Strength"))
summary(pureCFA.fit)

GetPattern(pureCFA.fit)

GetPrettyPattern(pureCFA.fit)

RMSEA(pureCFA.fit)

CheckMod(pureCFA.fit)

fit2 <- UseMod(pureCFA.fit)

RMSEA(fit2)

CheckMod(fit2)

fit3 <- UseMod(fit2)

RMSEA(fit3)

CheckMod(fit3)

GetPrettyPattern(fit3)

####

fit <- QuickEFAtoCFA(AD.R,3,1000)
GetPrettyPattern(fit)

fit <- QuickEFAtoCFA(AD.R,3,1000, cov.matrix=FALSE)
GetPrettyPattern(fit)

##########
R <- as.matrix(Harman74.cor$cov)
fit <- QuickJoreskog(R,4,145)
summary(fit)
GetPrettyPattern(fit)